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EXECUTIVE  SUMMARY 

Soldiers  of  the  1991  Gulf  War  were  given  the  drug  pyridostigmine  bromide  (PB)  for  pretreatment 
against  nerve  agents.  This  drug  has  since  been  associated  with  increased  incidence  and 
symptom  severity  of  Gulf  War  Illness.  The  goal  of  our  research  program  was  to  use  zebrafish 
as  a  model  organism  for  the  characterization  of  the  effects  of  pyridostigmine  bromide  on  gene 
expression  using  unbiased,  high-throughput  techniques,  specifically  mRNA-seq.  Additionally, 
we  intended  to  study  the  PB  response  in  the  presence  of  factors  that  may  modulate  the  effects 
of  PB  as  indicated  by  epidemiological  and  experimental  evidence,  such  as  genetic  variation  or 
stress.  Over  the  duration  of  the  award,  we  have  successfully  characterized  phenotypes 
associated  with  PB  exposure  and  measured  the  molecular  response  using  mRNA-seq.  We 
further  studied  the  effects  of  PB  in  the  context  of  three  different  genetic  backgrounds  and  under 
stressful  conditions.  Larval  zebrafish  were  utilized  to  study  the  systemic  response,  whereas 
results  from  experiments  with  adult  fish  provided  valuable  knowledge  of  effects  that  were  tissue- 
specific  and  those  that  were  observed  8  weeks  after  a  short-term  exposure.  Two  strains  of 
acetylcholinesterase  (AChE)  mutant  zebrafish  were  used  as  a  genetic  model  of  lowered  levels 
of  AChE  activity  and  two  other  AChE  inhibitors  were  studied  for  comparison  to  the  effects  to 
AChE  inhibition  by  PB.  These  experiments  laid  the  foundation  for  studying  the  off-target  effects 
of  PB  and  for  consideration  of  therapeutic  alternatives  to  PB,  respectively.  With  the  completion 
of  this  research,  we  have  provided  insight  into  the  predominant  changes  in  gene  expression  in 
response  to  PB  and  related  these  changes  to  the  phenotypes  of  PB-exposed  zebrafish.  We 
anticipate  that  a  greater  understanding  of  the  effects  of  PB  at  the  molecular  level  will  inform 
decisions  regarding  the  safety  and  use  of  PB  and  aid  in  the  development  of  therapies  for 
veterans  afflicted  by  Gulf  War  Illness. 


Contract  No.  W81XWH-09-1-0715  /  Proposal  #08273008 
HudsonAlpha  Institute  for  Biotechnology 
Richard  M.  Myers,  PI 
Final  report:  21  Sept  2009-20  Sep  2013 


Progress  and  expected  milestones: 


Build  small 
aquaculture  facility 


Technical  objective  2: 
Expose  fish  to  chemical 
and  screen  for  resistance 


Technical  objective  4: 
Assess  off-target  effects  by 
genetic  and  chemical  means 


Establish  gene  expression 
profiles  of  chemical  treatment 

with  RNA-seq 

100%  completed 


Assess  natural  variation  in 
response  to  chemical  and 
determine  effect  of  stress 
100%  completed 


4 


Award  number:  W81XWH-09-1-0715 

Title:  Applying  Genomic  and  Genetic  Tools  to  Understand  and  Mitigate  Damage  from  Exposure  to  Toxins 
Principal  Investigator:  Richard  M.  Myers,  Ph.D. 

INTRODUCTION 

Pyridostigmine  bromide  (PB)  was  given  to  approximately  250,000  soldiers  of  the  1991  Gulf  War 
as  a  prophylactic  against  potential  nerve  agent  attack  (7).  Pyridostigmine  is  a  peripheral 
acetylcholinesterase  (AChE)  inhibitor  that  reversibly  binds  AChE,  thereby  protecting  the  enzyme 
from  being  bound  by  nerve  agents,  which  are  irreversible  AChE  inhibitors.  Use  of  PB  has  since 
been  positively  associated  with  an  increased  incidence  of  Gulf  War  Illness  and  increased 
severity  of  the  symptoms  of  Gulf  War  Illness,  such  as  fatigue,  musculoskeletal  pain,  mood  and 
cognitive  disorders,  dermatological  irritation,  and  gastrointestinal  effects  (2).  Recent 
epidemiological  evidence  suggests  that  for  soldiers  in  the  most  forward  deployments,  PB 
exposure  was  the  greatest  factor  associated  with  Gulf  War  Illness  (3).  Research  also  suggests 
that  there  are  potential  interactions  between  the  effects  of  PB  and  additional  factors,  specifically 
stress  and  genetics.  Although  the  precise  mechanism  is  unclear,  treatment  with  PB  concurrent 
with  stressful  conditions  has  been  shown  to  result  in  aberrant  behavioral  patterns  and  altered 
gene  expression,  enzyme  activities,  and/or  reactive  oxygen  species  levels  in  the  brain  and 
muscle  (4-6).  Genetic  variation  has  also  been  associated  with  poor  health  in  Gulf  War  veterans. 
For  example,  variants  of  paraoxonase  1  and  butyrylcholinesterase  that  result  in  slower 
detoxification  of  organophosphates  and  PB  have  been  found  to  be  associated  with  an  increased 
likelihood  of  Gulf  War  Illness  symptoms  in  veterans  (7,  8).  Therefore,  it  was  our  objective  to 
study  the  effects  of  PB  alone  and  in  the  context  of  stress  and  varied  genetic  backgrounds  by 
using  zebrafish  as  a  model  organism  and  utilizing  the  high-throughput  technology  of  mRNA-seq. 


BODY 

The  goal  of  our  research  program  (W81XWH-09-1-0715)  was  to  utilize  genetic  and  genomic 
tools  to  determine  the  effects  of  PB  in  order  to  further  understand  the  pathogenesis  of  Gulf  War 
Illness.  Over  the  duration  of  the  award  we  completed  construction  of  an  aquaculture  facility  and 
advanced  knowledge  related  to  the  following  objectives:  i)  establishment  of  gene  expression 
profiles  of  chemical  treatment  with  RNA-seq,  ii)  determination  of  the  influence  of  natural 
variation  and  stress  on  the  response  to  chemical,  iii)  assessment  of  off-target  effects  by  genetic 
and  chemical  means,  and  iv)  exposure  of  fish  to  chemical  with  screening  for  resistance. 


Construction  of  an  aquaculture  facility 

Construction  of  the  aquaculture  facility  at  HudsonAlpha  necessitated  remodeling  of  a  room 
including  refinishing  of  the  floor  and  reconfiguration  of  drains  and  electrical  work.  Supplies  were 
ordered  and  the  aquaculture  facility  was  completed  with  the  installation  of  the  fish  racks  and 
water  system  taking  place  June  28-29,  2010.  The  final  inspection  by  the  HudsonAlpha  IACUC 
took  place  on  July  1,  2010.  All  personnel  that  have  had  direct  contact  with  the  animals 
completed  the  Working  with  Laboratory  Zebrafish  and  Working  with  the  IACUC  training  modules 
provided  by  the  American  Association  for  Laboratory  Animal  Science.  Following  inspection  and 
training,  final  approval  for  animal  protocols  was  received  from  the  HudsonAlpha  IACUC  and 
ACURO  on  July  9  and  July  26,  2010  respectively.  The  first  fish  were  received  from  the 
Zebrafish  International  Resource  Center  on  August  12,  2010.  The  facility  houses  5  strains  of 
zebrafish  including  the  AB  and  TU  strains,  which  are  robust  and  commonly-used  wild  type  fish, 
and  the  SJD  strain,  which  is  not  fully  inbred,  but  tends  to  have  less  genetic  variability  and  may 
be  more  amenable  for  genetic  and  genomic  experimentation.  We  also  acquired  the  achetm2a 
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and  achetm205  mutant  strains,  which  were  utilized  as  a  reference  for  the  phenotype  associated 
with  lowered  levels  of  AChE  activity. 


Technical  objective  1:  establishment  of  gene  expression  profiles  of  chemical  treatment 
with  RNA-seq 

Dose-response  and  time-course  experiments  in  larval  and  adult  zebrafish  were  conducted  to 
determine  alterations  in  gene  expression  profiles  in  response  to  PB  exposure.  In  addition  to 
generating  preliminary  data  on  the  effects  of  PB,  these  studies  were  essential  for  identifying  the 
optimal  dosage  and  duration  of  treatment  for  subsequent  experiments. 

During  the  Gulf  War,  PB  was  administered  in  blister  packs  of  21  pills,  with  each  pill  containing  a 
30  mg  dose  of  PB  to  be  taken  every  eight  hours  while  the  threat  of  nerve  attack  persisted  ( 1 ). 
This  dosage  is  expected  to  maintain  sustained  inhibition  of  AChE  activity  by  approximately  30%. 
To  examine  the  dose-response  to  PB  in  larval  zebrafish  and  identify  the  proper  dose  for 
subsequent  experiments,  we  conducted  a  study  in  which  pooled  TU  zebrafish  larvae  (n  =  20) 
were  treated  with  solutions  of  0,  0.01,  0.1,  1,  or  10  mmol  PB/L  (n  =  4  per  treatment)  from  4-7 
days  post-fertilization  (dpf),  the  period  immediately  following  completion  of  organogenesis  and 
hatching  of  the  embryo.  Following  treatment,  the  larvae  were  rinsed  three  times  with  ice-cold 
phosphate-buffered  saline  (PBS)  and  snap  frozen  on  dry  ice.  Ten  larvae  were  preserved  in 
Trizol  for  extraction  of  nucleic  acids  for  mRNA-seq,  while  the  remaining  larvae  were  used  for 
preparation  of  homogenates  for  enzymatically-active  protein  using  the  PARIS  Kit  (Ambion,  Inc.). 
Phenotypes  associated  with  the  primary  action  of  PB,  altered  acetylcholine  signaling,  were 
characterized  by  determination  of  AChE  activity  levels  and  assessment  of  body  posture 
following  movement  elicited  by  the  startle  response.  AChE  activity  was  determined  using  the 
AmplexRed  Acetylcholinesterase  Assay  (Invitrogen,  Inc.).  Individual  zebrafish  larvae  were 
treated  with  0,  0.01,  or  10  mmol  PB/L  (n  =  4  per  treatment)  for  72-h  and  movement  was  video 
recorded  during  trials  in  which  the  startle  response  was  elicited  ten  times  by  approaching  the 
larvae  with  a  blunt  object.  The  angle  of  the  tail  bend  at  rest  following  the  startle  response  was 
measured  as  an  indicator  of  the  neuromuscular  effects  of  PB. 


Inhibition  of  AChE  activity  was  dose- 
dependent  in  larval  zebrafish,  with 
doses  of  0.1  and  1  mmol  PB/L 
producing  levels  of  AChE  inhibition  that 
were  within  the  therapeutically  relevant 
range  of  -20-50%  (Figure  1).  We  also 
observed  neuromuscular  effects  of  PB 
in  fish  treated  with  10  mmol  PB/L  as 
evidenced  by  a  greater  angle  in  the  tail 
bend  at  rest  following  the  startle 
response  in  zebrafish  larvae  (Figure  2, 
p  =  0.01).  The  increased  bend  of  the 
tail  suggests  increased  contractile 
tension,  which  may  be  due  to 
hyperstimulation  of  acetylcholine 
receptors  in  response  to  inhibition  of 
the  AChE  activity  and  the  resulting 
accumulation  of  acetylcholine  at  the 


Figure  1.  Inhibition  of  acetylchol  nesterase  by  pyrdostigmme  bromide  is  dose- 
dependent  in  larval  zebrafish.  Data  a-e  means  (n  =  4)  ±  s.e. 
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neuromuscular  synapse.  These  data  indicate  that  zebrafish  larvae  are  an  appropriate  model  for 
the  study  of  PB:  we  were  able  to  achieve  levels  of  inhibition  similar  to  that  used  therapeutically; 
and  we  observed  neuromuscular  effects  at  a  higher  dose,  which  is  in  line  with  the  results  that 
have  been  reported  previously  in  other  species  (7). 


PyrkJosDgmine  brande  (mmoVL) 

Figure  2.  72-hr  treatment  of  TU  larval  zebrafish  with  pyridostigmine  bromide  produces  subtle  phenotypic  effects  after 
the  startle  response. 


Preliminary  mRNA-seq  data  from  two  samples  per  treatment  dose  was  produced  to  assess  the 
effects  of  PB  at  the  molecular  level.  RNA  isolated  using  the  Trizol  protocol  was  used  as  input 
for  the  generation  of  libraries  via  a  protocol  developed  by  our  lab  that  uses  Nextera  technology 
(9).  The  libraries  were  sequenced  on  an  lllumina  HiSeq  and  the  resulting  sequence  was  aligned 
to  zebrafish  RefSeq  exons  (zv9)  using  Bowtie  (70).  Variance  stabilization  of  the  data  and 
differential  expression  for  each  PB-treated  group  compared  to  controls  were  calculated  using 
statistical  methodologies  developed  in  our  lab,  as  described  by  Reddy  et  al  (7  7).  This  analysis 
identified  81  transcripts  differentially  expressed  in  response  to  10  mmol/L  PB  treatment  when 
compared  to  control  larvae;  there  were  four  and  three  differentially  expressed  transcripts  in 
response  to  the  1  and  0.1  mmol/L  PB  treatments,  respectively.  The  DAVID  Functional 
Annotation  Tool  (72,  73)  was  used  to  test  for  functional  annotation  enrichments  within  the 
differentially  expressed  transcripts  at  10  mmol  PB/L,  with  significant  enrichment  assigned  to 
clusters  with  scores  greater  than  or  equal  to  1.30  or  terms  meeting  an  FDR  cutoff  of  5.  The  44 
upregulated  transcripts  were  significantly  enriched  for  hemopexin/matrixin  repeat,  peptidase 
activity,  and  endopeptidase  inhibitor  activity,  with  data  also  suggestive  of  enrichment  for  heme 
binding  (FDR  =  7.15).  The  37  downregulated  transcripts  were  significantly  enriched  for 
cognition  and  sensory  perception,  intermediate  filament  proteins,  EF-HAND  and  calcium 
binding,  with  data  suggestive  of  enrichment  for  the  KEGG  pathway  related  to  biosynthesis  of 
unsaturated  fatty  acids  (FDR  =  8.98). 

The  order-restricted  inference  for  ordered  gene  expression  (ORIOGEN)  software  (14)  was  used 
to  look  for  dose-responsive  genes.  Although  we  were  unable  to  obtain  statistical  significance  at 
the  genome-wide  level  with  this  small  preliminary  data  set,  enrichment  of  the  genes  for  which  a 
significant  dose-response  was  found  at  the  gene  level  (p  <  0.05)  identified  terms  that  appear 
plausible  based  on  the  known  molecular  action  of  PB.  At  the  gene-level,  946  genes  were  found 
to  have  an  increasing  profile  and  597  genes  had  a  decreasing  profile.  Decreasing  transcripts 
were  enriched  for  terms  related  to  transmission  of  nerve  impulses,  neuron  projection 
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development,  calcium  channels  and  binding  proteins,  intermediate  filament,  immunoglobulin, 
and  metabolism.  Increased  transcripts  were  enriched  for  peptidase  activity,  GTP  binding, 
transmembrane  proteins,  cell  cycle,  interferon  regulatory  factor,  protein  localization,  FAD- 
dependent  pyridine  nucleotide-disulphide  oxidoreductase,  and  endoplasmic  reticulum.  The 
results  of  both  analyses  are  indicative  of  PB  having  molecular  effects  related  to  its  primary 
mechanism  of  action,  as  evidenced  by  enrichment  for  terms  associated  with  the  synapse  and 
regulation  of  calcium  levels.  However,  there  was  also  enrichment  for  terms  related  to  structural 
components  of  the  cell,  the  immune  response,  and  peptidase  activity  that  may  also  be 
informative. 

The  results  of  the  dose  response  study  provided  valuable  information  regarding  the  appropriate 
dosage  to  mimic  therapeutic  conditions.  There  is  epidemiological  evidence  that  consumption  of 
an  increased  number  of  PB  pills  was  associated  with  increased  risk  of  Gulf  War  Illness  and 
severity  of  symptoms  (2,  15).  We  designed  a  study  to  examine  the  effects  of  PB  treatment  after 
exposure  times  ranging  from  4  hours  to  7  days,  based  on  a  report  that  of  the  soldiers  that 
recalled  using  PB,  the  majority  took  PB  for  less  than  seven  days  (76).  Pooled  TU  zebrafish 
larvae  (n  =  20)  were  exposed  to  0.5  mmol/L  PB  solution  for  4,  24,  72,  or  168-h  beginning  at  4 
dpf;  system  water  was  used  as  a  control.  Sample  preparation,  AChE  activity  assessment,  and 
mRNA-seq  were  performed  as  described  for  the  dose-response  study. 

Our  analysis  of  the  effects  on  AChE  activity 
took  into  account  the  exposure  to  PB,  the 
duration  of  exposure  and  a  possible 
interaction  between  PB  and  exposure  time. 

Pyridostigmine  significantly  inhibited  AChE 
activity  (p  <  0.01),  however  there  was  no 
effect  of  time  or  an  interaction  between  PB 
and  the  duration  of  exposure  (Figure  3). 

Three  samples  per  treatment  were 
prepared  for  mRNA-seq;  sequencing 
generated  -18-28  million  reads  aligning  to 
exons  of  the  zebrafish  transcriptome. 

Analysis  using  a  linear  regression  model 
with  terms  for  PB  treatment,  time  and  the 
treatment  x  time  interaction  revealed  a 
significant  effect  of  time,  with  7,322 
transcripts  identified  as  differentially 
expressed.  These  transcripts  were 
enriched  for  broad  categories  such  as  DNA  metabolism,  homeostatic  process  and  cell  cycle.  In 
this  model,  there  was  not  a  significant  effect  of  PB  treatment  nor  a  treatment  x  time  interaction 
for  any  transcripts.  Data  for  each  time  point  was  analyzed  separately  (control  vs.  PB-treated), 
identifying  136,  327,  980,  and  460  DE  genes  at  4,  24,  72,  and  168  hours  respectively.  At  four 
hours,  there  was  significant  enrichment  for  DNA-dependent  DNA  replication  and  DNA  metabolic 
processes.  At  24  hours,  there  was  significant  enrichment  for  response  to  metal  ion  or  inorganic 
substance,  basic-leucine  zipper  transcription  factor,  and  WD40  repeat.  At  72  hours,  there  was 
significant  enrichment  for  cytoskeleton,  proteasome  complex,  methylation,  and  endoplasmic 
reticulum.  At  168  hours,  there  was  significant  enrichment  for  establishment  of  protein 
localization,  cell  redox  homeostasis,  intracellular  signaling  cascade,  and  endoplasmic  reticulum. 


Figure  3.  Acetylcholinesterase  s  significantly  inh  bited  by  PB  exposure, 
whereas  there  is  no  effect  of  time.  Data  are  means  (n  =  3)  ±  s.e. 
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The  dose  response  and  time  course  studies  provided  valuable  data  regarding  the  response  of 
larval  zebrafish  to  PB  at  the  phenotypic  and  molecular  levels.  Previous  studies  have  suggested 
that  the  response  to  pharmaceuticals  or  toxins  is  very  similar  in  larvae  and  adults  when  using 
zebrafish  as  a  model  organism  (77).  However,  to  verify  this  result  for  PB,  we  conducted  a 
preliminary  study  of  the  response  to  PB  in  adult  TU  zebrafish.  Adult  zebrafish  (8  months  of  age) 
were  removed  from  the  communal  tank  and  placed  into  individual  beakers  with  300  ml.  of  fish 
water  and  allowed  to  acclimate  for  3  days  before  the  initiation  of  a  24-h  treatment  with  solutions 
of  0,  0.1,  1,  10,  100,  or  1000  pmol  PB/L.  One  male  and  one  female  fish  were  treated  at  each 
dose.  At  the  end  of  the  treatment  period,  the  fish  were  euthanized  by  submersion  in  ice  water 
followed  by  decapitation  with  a  razor  blade.  The  head  and  the  body  of  each  fish  were 
separately  prepared  for  isolation  of  protein  and  RNA  for  the  measurement  of  AChE  activity  and 
mRNA-seq  by  the  methods  as  described  for  the  larval  studies.  Alignment  to  the  RefSeq  exons 
of  the  zebrafish  genome  resulted  in  21-48  million  reads  per  sample  for  downstream  analyses. 

Consistent  with  reports 
that  PB  is  unable  to 
cross  the  blood-brain 
barrier  (78),  there  were 
no  significant  effects  of 
PB  treatment  on  AChE 
activity  in  the  head. 

However,  AChE  activity 
in  the  body  was 
significantly  inhibited  in 
adult  zebrafish  treated 
with  either  10  or  1000 
pmol  PB/L  relative  to 
controls  (Figure  4). 

Taken  together  with 
observations  that  the 
dose-response  to  AChE 
activity  was  much 
stronger  in  the  male  fish 
than  the  female  fish  (R2 
=  0.71  vs.  0.14),  we  chose  to  do  mRNA  profiling  by  mRNA-seq  using  samples  from  the  heads 
and  bodies  of  the  male  fish  treated  with  0,  10,  and  1000  pmol  PB/L.  Using  unsupervised 
hierarchical  clustering  based  on  the  293  most  variable  genes,  samples  clustered  first  by  the 
region  of  the  body  (Figure  5).  Within  the  tree  node  representing  the  body  samples,  the  PB- 
treated  samples  clustered  away  from  the  control,  whereas  the  differences  were  much  smaller 
among  the  samples  from  the  head,  with  no  separation  of  the  PB-treated  samples.  Our  statistical 
analysis  of  all  expressed  genes  confirmed  that  the  largest  effects  were  due  to  the  region  from 
with  the  sample  originated,  with  2,884  transcripts  being  differentially  expressed  between  the 
head  and  the  body  (FDR  <  0.05).  As  expected,  transcripts  that  were  significantly  differentially 
expressed  with  at  least  50%  greater  expression  in  the  head  were  enriched  for  genes  associated 
with  the  eyes  and  brain,  whereas  the  enriched  ZFIN  anatomy  terms  associated  with  greater 
expression  in  the  body  included  trunk  musculature  and  the  whole  organism.  Likely  due  to 
current  lack  of  replicates,  there  were  no  transcripts  that  met  this  strict  criteria  (FDR  <  0.05)  for 
the  main  effect  of  PB  treatment  or  for  the  treatmenhtissue  interaction.  However,  407  mRNAs 
were  affected  by  PB  treatment  and  1086  mRNAs  had  a  response  to  PB  that  also  depended  on 
the  origin  of  the  sample  at  the  gene-level  (p  <  0.05).  Enrichment  of  terms  for  the  PB  response 
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Figure  4.  Acetylcholinesterase  activity  in  the  head  or  body  of  adult  zebrafish  treated  with 
pyridostigmine  bromide  for  24-h.  Data  are  means  t  s.e. 
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included  calcium  ion  binding,  regulation  of  peptidase  activity,  positive  regulation  of 
macromolecular  biosynthetic  process,  and  sphingolipid  metabolism.  Enrichment  of  terms  for  the 
PB  by  tissue  interaction  included  DNA  metabolic  processes,  cell  cycle,  cytoskeleton,  cell 
projection  part,  sensory  organ  development,  DNA  replication,  DNA-dependent  ATPase  activity, 
tetrtricopeptide  region,  DNA  binding,  and  IQ-calmodulin  binding  region.  In  an  analysis  of  the 
samples  from  the  head  or  the  body  for  the  effects  of  PB  separately,  there  were  893  transcripts 
differentially  regulated  by  PB  in  the  body  and  639  transcripts  identified  as  differentially 
expressed  in  the  head.  Terms  enriched  in  body  DE  genes  included  transcription  factor  activity, 
microtubule  cytoskeleton,  sphingolipid  metabolism.  Terms  enriched  in  the  head  DE  genes 
included  transcription  factor  complex,  ribosome,  endosome,  intracellular  non-membrane- 
bounded  organelle,  positive  regulation  of  gene  expression  (protein  heterodimerization  activity 
and  nuclear  translocation),  sensory  organ  development,  ferritin,  and  chromosome  localization. 

These  studies  laid  the  foundation  for  the  more  complex  studies  of  the  PB  response  within  the 
context  of  other  factors  in  both  larval  and  adult  zebrafish.  In  addition  to  establishment  of  the 
appropriate  dose  and  duration  for  the  PB  treatments,  these  preliminary  mRNA-seq  data  sets 
were  suggestive  of  PB  effects  at  genes  related  to  the  synapse  and  regulation  of  calcium  levels, 
as  well  as  additional  targets  related  to  cytoskeletal  structures,  regulation  of  protein  digestion  and 
degradation,  and  regulation  of  cellular  homeostasis,  particularly  redox  and  iron  ion  homeostasis. 
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Figure  5.  Unsupervised  hierarchical  cluster  of  RNA-seq  data  from  the  head  and  bodies  of  adult  zebrafish 
treated  with  pyridostigmine  bromide. 


Technical  objective  2:  determination  of  the  influence  of  natural  variation  and  stress  on 
the  response  to  chemical 

Effects  of  genetic  variation  on  the  response  to  pyridostigmine  bromide 

We  conducted  an  experiment  with  the  goal  of  further  characterizing  the  dose-dependent  effects 
of  PB  as  well  as  determining  the  effect  on  natural  genetic  variation  on  the  response  to  the  drug. 
Pooled  groups  of  20  larvae  from  each  of  three  common  laboratory  strains  -  AB,  TU,  and  SJD  - 
were  treated  with  solutions  of  0,  0.001,  0.01,  0.1,  1,  or  10  mmol  PB/L  for  72-h.  The  treatment 
period  corresponded  to  the  period  of  4-7  days  post-fertilization.  At  7  dpf,  larvae  were 
euthanized  and  samples  were  prepared  as  described  for  previous  larval  studies. 
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Phenotypes  associated  with  PB  treatment  were  assessed  by  determination  of 
acetylcholinesterase  activity  and  by  using  motility  and  neurobehavioral  assays. 

Acetylcholinesterase  (AChE)  activity  was 
measured  by  the  Amplex  Red 
Acetylcholinesterase  Assay  (Invitrogen) 
for  determination  of  the  efficacy  of 
treatment  and  to  identify  the  clinically 
relevant  doses.  Pyridostigmine  bromide 
treatment  significantly  inhibited  AChE 
activity  in  a  dose-dependent  manner  in 
all  three  strains  of  zebrafish.  There  was 
no  significant  effect  of  strain  nor  a 
treatment  by  strain  interaction.  The  dose- 
dependent  effect  of  PB  exposure 
produced  inhibition  at  levels  significantly 
different  from  the  control  at  all  doses 
greater  than  or  equal  to  0.01  M/L  (Figure 
6).  For  the  motility  analysis  zebrafish 
larvae  were  placed  in  24-well  plates  (1 
fish  per  well)  and  the  PB  and  control 
solution  treatments  were  randomized 
across  each  plate  (Figure  7).  A  total  of  four  plates  per  strain  were  used  for  this  experimental 
series,  with  two  plates  receiving  the  stress  treatment  and  two  plates  serving  as  controls. 
Patterns  of  movement  were  recorded  and  quantified  using  the  ZebraLab  equipment  and 
VideoTrack  software  (Noldus  Technologies).  The  protocol  tracked  movement  for  five  minutes 
under  normal  light  conditions,  an  abrupt  transition  to  five  minutes  of  darkness,  and  an  additional 
five  minutes  under  normal  light  conditions  post-stimulus.  Movement  patterns  in  the  whole  well 
and  the  center  of  the  well  (50%  of  total  area)  were  used  for  determining  patterns  in  total 
locomotion  as  well  as  place  preference,  both  of  which  can  be  used  as  a  behavioral  indicator  of 
anxiety  (19).  This  assay  is  similar  in  concept  to  the  open  field  test  used  in  rodents  and  is  very 
similar  to  a  protocol  that  has  recently  been  published  by  Schnorr  et  al.  (20),  although  our 
protocol  also  includes  characterization  of  post-stimulus  behavior. 
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Figure  6.  Pyridostigmine  bromide  Inhibits  acetylcholinesterase  activity  in  three  strains 
of  zebrafish  in  a  dose-dependent  manner.  Data  are  means  (n  =  3)  ±  s.e. 
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Figure  7.  Motility  and  behavioral  analysis  of  zebrafish  larvae  under  baseline  conditions  and  the 
response  to  abruptly  changing  light  conditions. 
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Figure  8.  Strain-specific  c  defences  in  mobiity  and  place  preference  dunng  {A)  base  ne  light  phase, 

{8)  dark  stimulus,  and  (C)  post-stimulus  light  phase.  Data  are  means  (n  =  47-48}  ±  s.e. 

There  were  a  number  of  differences  in  locomotion  and  place  preference  noted  between  strains, 
with  differences  occurring  during  all  phases  of  the  assay  (Figure  8).  PB  treatment  resulted  in 
decreased  distance  traveled  in  the  center  at  baseline  and  during  the  dark  phase,  decreased 
total  duration  in  center  during  dark  (trend  at  baseline),  decreased  percentage  of  distance 
traveled  in  center  during  dark  (trend  at  baseline),  as  well  as  a  trend  towards  fewer  entries  into 
center  at  baseline  (Figure  9).  A  treatment  by  strain  interaction  was  observed  for  total  distance 
at  baseline,  distance  and  percent  total  distance  in  center  after  stimulus.  Analysis  of  each  strain 
individually  for  strain-specific  effects  of  PB  identified  the  center  distance  and  duration  (alone  and 
as  percent  of  total)  significantly  decreased  in  dark  for  the  AB  strain,  no  differences  for  the  SJD 
strain,  and  in  the  TU  strain,  there  were  trends  for  decreased  total  distance,  duration,  and 
average  velocity  in  the  center  at  baseline,  as  well  as  a  decreased  number  of  center  entries  in 
dark.  In  concordance  with  data  from  another  study  utilizing  the  light:dark  test  in  assessment  of 
chemical  response  (ethanol)  in  different  strains  of  zebrafish  (21),  we  also  found  that  there  were 
significant  differences  in  the  movement  patterns  of  the  three  strains  of  zebrafish,  but  despite 
these  differences,  we  were  able  to  detect  a  response  to  PB  treatment.  The  distance  and 
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Figure  9.  PB  effects  on  mobMy  and  behavior  by  stran  during  (A)  base^e  light  phase. 
(8)  dark  stimulus,  and  (C)  post-sbrnulus  light  phase.  Data  are  means  (n  =  23-24)  t  s.e. 


duration  of  movement  in  the  center  of  the  well  was  decreased  by  exposure  to  PB,  suggesting 
that  PB  has  anxiogenic  effects  in  zebrafish  larvae. 

For  the  molecular  characterization  of  the  PB  response  in  three  strains  of  zebrafish,  54  samples 
(n=3  per  treatment  per  strain)  were  used  for  mRNA-seq.  RNA  samples  were  prepared  into 
libraries  for  mRNA-seq  using  the  Nextera  protocol  (9)  and  sequenced  as  paired-end  50-bp  on 
an  lllumina  HiSeq.  With  upgrades  to  our  analysis  pipeline  following  our  preliminary 
experiments,  we  completed  our  alignments  to  the  Ensembl  database  using  TopHat  (22)  and 
compared  the  aligned  transcriptome  to  the  reference  genome  and  RNAs  present  in  the  Refseq, 
Ensembl,  and  UCSC  mRNA  databases  for  identification  of  changes  in  the  transcript  sequence  - 
such  as  indels,  single-  and  multiple-nucleotide  changes  using  snpEff  (23).  From  the  alignments, 
HTseq  was  utilized  to  generate  raw  counts  for  Ensembl  genes  (ens67).  The  analysis  and 
interpretation  of  RNA-seq  data  may  vary  by  the  methodology  used,  so  we  compared  several  of 
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the  most  commonly  used  statistical  packages  for  determination  of  differential  expression  such 
as  DESeq,  edgeR,  and  SAMseq  (24-26).  For  the  questions  of  interest  in  these  experiments,  our 
analysis  indicated  that  DESeq  was  the  most  appropriate  tool  because  SAMseq  cannot  handle 
multivariate  models  or  interactions  and  edgeR  had  greater  variability  in  the  estimates  of 
dispersion  than  DESeq.  Therefore,  DESeq  was  used  for  variance  stabilization  and  differential 
expression  calculations  for  the  count  data  from  HTSeq.  With  a  threshold  for  the  false  discovery 
rate  of  0.05,  we  found  356,  10,416,  and  87  changes  in  gene  expression  attributable  to  the 
effects  of  PB,  strain  and  PB  x  strain  interaction,  respectively.  We  used  the  DAVID  Functional 
Annotation  Tool  (12,  13)  to  characterize  the  enrichment  of  altered  transcripts.  This  identified  the 
terms  associated  with  PB  treatment  as  endopeptidase  activity,  aminoacyl-tRNA  biosynthesis, 
actin,  and  negative  regulation  of  apoptosis  (p53  signaling  pathway).  There  was  significant 
enrichment  for  terms  associated  with  iron  ion  binding,  DNA-dependent  regulation  of 
transcription,  ribonucleotide  binding,  and  alpha-linolenic  acid  metabolism  in  the  genes  altered  by 
the  PB  x  strain  interaction.  The  dose-response  to  PB  was  also  assessed  using  the  ORIOGEN 
tool  with  the  variance  stabilized  data  produced  by  DESeq.  For  this  analysis,  samples  within 
each  strain  were  binned  based  on  the  AChE  activity  for  the  sample:  control,  PB-treated  with  less 
than  20%  inhibition,  20-50%  inhibition,  and  greater  than  50%  AChE  inhibition.  Using  ORIOGEN 
to  study  the  dose-response  using  this  grouping  identified  24  transcripts  with  a  decreasing  profile 
and  16  transcripts  that  were  upregulated  with  increasing  degrees  of  inhibition  of  AChE.  The 
DAVID  Functional  Analysis  Tool  is  generally  not  appropriate  for  gene  lists  this  small;  however, 
the  Integrative  Multi-species  Prediction  (IMP)  tool  (27,  28)  was  developed  specifically  for  gene 
sets  of  this  size  and  was  used  to  search  for  enrichments  of  biological  processes  predicted  to  be 
involved  with  these  genes.  The  upregulated  genes  were  enriched  for  antigen  processing  and  t- 
cell  mediated  cytotoxicity,  whereas  the  downregulated  genes  were  enriched  for  neuron  synaptic 
transmission,  positive  regulation  of  G  protein  coupled  receptor  protein  signaling  pathway, 
circadian  sleep  wake  cycle,  regulation  of  renal  sodium  excretion,  and  elevation  of  cytosolic 
calcium  ion  concentration.  Enrichment  of  biological  processes  for  genes  most  downregulated  at 
or  below  therapeutic  levels  of  AChE  inhibition  also  included  dendrite  guidance,  visceral  motor 
neuron  differentiation,  and  muscle  cell  fate;  upregulated  genes  were  enriched  for  carbohydrate 
transport,  peptidyl  serine  phosphorylation,  long  term  memory,  microtubule  anchoring,  long  term 
potentiation,  antibacterial  humoral  response,  and  regulation  of  microtubule  polymerization. 

Due  to  the  large  number  of  genes  identified  as  statistically  different  between  strains,  subsequent 
analyses  were  undertaken  within  each  strain  to  identify  dose-dependent  effects  and  those 
effects  specifically  associated  with  treatments  near  the  clinically  relevant  range  of  -20-50% 
inhibition  of  AChE  activity  (7).  Similar  to  the  results  for  the  motility  assay,  there  was  the  least 
significant  response  in  the  SJD  strain,  whereas  analysis  of  the  TU  and  AB  strains  produced 
larger  sets  of  differentially  expressed  genes  (Figure  10).  Expression  of  87  genes  was  altered  in 
the  SJD  strain,  with  enrichment  for  peptidases.  Exposure  to  PB  altered  expression  of  235 
transcripts  in  AB  larval  zebrafish,  with  enrichment  for  dynamin,  organic  acid  biosynthetic 
process,  oxidation-reduction,  proteasome,  and  interferon  regulatory  factor.  In  the  TU  strain 
there  was  enrichment  for  response  lectin/sugar  binding,  apoptosis,  and  response  to  bacterium 
in  the  316  differentially  regulated  transcripts.  The  greatest  degree  of  overlap  in  differentially 
expressed  genes  was  found  between  the  AB  and  TU  strains  (Figure  10)  and  within  those  70 
genes,  there  was  overrepresentation  within  the  interferon  regulatory  factors  and  the  p53 
signaling  pathway.  The  mRNA-seq  data  from  this  experiment  confirmed  the  findings  of  the 
preliminary  studies  with  enrichments  for  both  neurological  function  and  calcium  signaling,  as 
well  as  for  cytoskeletal  structures,  regulation  of  protein  digestion  and  degradation,  and 
regulation  of  cellular  homeostasis,  particularly  redox  and  iron  ion  homeostasis. 
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Figure  10.  Overlap  in  genes  differentially  expressed  in  response  to  PB  in  Figure  11.  Overlap  in  SNPs  identified  by  snpEff  from  RNA-seq  data 

larval  zebrafish  of  three  strains.  from  three  strains  of  zebrafish. 

We  have  also  begun  exploring  the  data  provided  by  the  output  of  SnpEff  (21),  a  SNP  discovery 
and  annotation  tool.  Using  these  algorithms,  we  identified  14,018  nonsynonymous  nucleotide 
changes  that  might  alter  protein  function.  Of  these  SNPs,  the  fewest  SNPs  were  observed  in 
the  TU  strain  and  the  greatest  number  of  SNPs  was  found  in  the  SJD  strain  (Figure  11).  The 
relative  number  of  SNPs  observed  and  extent  of  overlap  between  strains  for  differential 
expression  and  SNPs  is  consistent  with  our  expectations  because  the  TU  strain  was  the  strain 
used  for  the  reference  genome  and  the  AB  and  TU  strains  are  more  genetically  similar  than  the 
SJD  strain  (29).  It  is  notable  that  there  is  little  overlap  in  the  PB-dose  responsive  genes 
between  strains  with  over  70%  of  PB-responsive  genes  being  unique  to  one  strain.  Although 
there  is  a  greater  degree  of  overlap  in  SNPs,  no  less  than  40%  of  SNPs  identified  by  strain  are 
unique  for  that  strain.  Taken  together  with  the  motility  data,  this  provides  an  excellent 
opportunity  to  understand  genotype-driven  strain-specific  differences  in  gene  expression, 
generally  and  specifically  in  response  to  PB  exposure.  The  Fisher’s  exact  test  found  that  there 
was  not  a  greater  than  expected  number  of  SNPs  within  the  differentially  expressed  genes 
within  any  of  the  strains.  However,  we  anticipate  that  we  will  use  the  Gene  Set  Association 
Analysis-SNP  (GSAA-SNP)  software  (30)  to  further  investigate  possible  relationships  between 
the  SNPs  identified  within  each  strain  and  the  differentially  expressed  gene  data  sets.  This 
work  could  be  translated  into  a  greater  ability  to  predict  whether  there  are  additional  human 
populations  that  might  be  more  susceptible  to  the  effects  of  PB  and  whether  the  response  to  PB 
varies  by  genotypic  differences  at  other  loci  beyond  those  for  detoxification  enzymes,  such  as 
butyrylcholinesterase  and  paraoxonase. 


Effects  of  stress  on  the  response  to  pyridostigmine  bromide 

In  addition  to  exposure  to  PB,  soldiers  of  the  Gulf  War  of  1991  were  also  exposed  to  a  wide 
variety  of  additional  stressors.  Given  evidence  for  crosstalk  between  the  acetylcholine  and 
stress  signaling  pathway  via  a  number  of  molecular  mechanisms  (6,  31-33),  we  sought  to 
identify  the  effects  of  PB,  with  or  without  concurrent  stress,  in  larval  zebrafish  of  the  AB  strain. 
We  conducted  a  study  with  the  aim  of  understanding  the  effects  of  PB  exposure  and  stress, 
alone  or  in  combination,  in  larval  zebrafish.  Pooled  groups  of  AB  zebrafish  larvae  (n=20  larvae) 
were  randomly  assigned  to  one  of  four  treatment  groups:  control,  PB,  stress,  or  PB+stress 
(n=24  per  group).  PB-treated  fish  were  kept  in  0.5  mmol  PB/L  for  72  hours  from  4-7  days  post- 
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Table  1.  The  effects  of  PB  on  acetylcholinestrase  activities  in  larval 
zebrafish  were  not  affected  by  stress.  Data  are  means  (n  =  5)  ±  s.e. 

Treatment 

Acetylcholinesterase  activity  (%  control) 

control 

100.0  ±3.0 

stress 

96.6  ±4.0 

PB 

56.0  ±  2.6  * 

PB+stress 

54.2  ±  3.2  * 

•indicates  significant  difference  from  control  (p<0.05) 
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Figure  12.  MoMrty  and  cxace  {deference  are  altered  by  PB  and  stress  in 
larval  zebraf  sh  during  (A)  Daselme  light  phase.  (B)  darn  stimulus,  and  (C) 
post-stimulus  light  phase.  Data  are  means  (n  =  31-36}  ±  s.e. 


fertilization.  To  induce  a 
stressful  condition,  larvae  were 
swirled  for  1  min  twice  per  day 
starting  at  4  dpf  (34).  Samples 
were  collected  for  protein  for 
enzymatic  assays  and  for 
nucleic  acid  extraction  at  7  dpf 
as  described  for  the  previous 
larval  experiments,  motility  and 
behavior  assays  were  also 
performed  as  described  in  Figure  7. 
Acetylcholinesterase  activity  was 
inhibited  to  within  the  clinically  relevant 
range  of  20-50%  by  PB  exposure, 
whereas  stressful  conditions  had  no 
effect  on  acetylcholinesterase  activity 
(Table  1). 

A  critical  component  of  our  research  into 
the  effect  of  PB  is  to  determine  whether 
the  drug  produces  a  phenotype  that  is 
comparable  to  the  symptoms  of  Gulf 
War  Illness.  In  this  experiment,  PB- 
treated  larvae  were  most  affected  post¬ 
stimulus,  with  a  lower  percent  of  total 
distance  traveled  in  the  center  and  a 
lower  duration  of  time  spent  in  the  center 
after  the  changing  of  light  conditions, 
suggesting  that  PB-treated  fish  were 
more  anxious  (Figure  12).  Conversely, 
during  the  dark  phase  stress-treated 
larvae  tended  to  travel  a  greater 
distance  in  total  and  also  a  greater 
proportion  of  that  distance  tended  to  be 
in  the  center  of  the  well,  suggesting  that 
inhibitions  were  decreased  in  stressed 
fish.  When  the  data  was  binned  by 
velocity,  there  were  more  high  velocity 
movements  in  alarm  pheromone- 
exposed  fish  during  the  dark  phase, 
suggesting  more  erratic  movement 
patterns  induced  by  stressful  conditions. 

Nucleic  acids  were  isolated  from  20 
samples  (n=5  per  group).  We 
sequenced  the  RNAseq  libraries  as  50- 
bp  paired  end  reads  on  an  lllumina 
HiSeq,  producing  66-73  million  reads 
aligning  to  the  zebrafish  transcriptome 
for  each  sample.  Using  our  analysis 


16 


Award  number:  W81XWH-09-1-0715 

Title:  Applying  Genomic  and  Genetic  Tools  to  Understand  and  Mitigate  Damage  from  Exposure  to  Toxins 
Principal  Investigator:  Richard  M.  Myers,  Ph.D. 


pipeline  as  described  for  the  natural  genetic  variation  study  we  identified  151  genes  differentially 
expressed  in  response  to  PB  treatment,  358  genes  differentially  expressed  under  stressful 
conditions,  and  57  genes  for  which  there  was  an  interaction  between  PB  treatment  and  stress. 
Functional  annotation  clustering  using  the  DAVID  analysis  tool  revealed  that  within  genes 
downregulated  by  PB  there  was  overrepresentation  of  annexin,  phospholipid  binding,  calcium 
binding,  EF-hand  motif,  and  actin,  and  within  genes  upregulated  by  PB  there  was  enrichment  for 
sequence-specific  DNA  binding.  Transcripts  altered  by  exposure  to  stress  were  enriched  for  the 
respiratory  chain,  ribosome,  response  to  wounding,  gelsolin  and  fibrinogen,  interferon  regulatory 
factor,  thick  filaments,  cell-cell  binding  and  calcium/phospholipids-binding  and  annexin.  An 
interaction  between  PB  and  stress  was  associated  with  genes  for  extracellular  regions, 
specifically  extracellular  matrix  proteins.  There  was  an  overlap  of  77  genes  that  were 
dysregulated  by  PB  and  stress,  with  no  interaction;  response  to  bacterium  and  calcium  signaling 
were  significantly  enriched  in  these  genes.  These  results  confirm  the  effects  of  PB  on  genes 
associated  with  calcium  signaling  and  cytoskeletal  structures.  Furthermore,  we  observed  an 
interaction  between  PB  and  stress  on  genes  for  extracellular  matrix  proteins,  which  are  also 
important  for  the  determination  of  cell  structure  and  function. 

The  studies  of  PB  and  stress  exposures  in  larval  zebrafish  were  useful  for  identifying  large 
systemic  changes  in  gene  expression  which  might  be  relevant  to  Gulf  War  Illness,  a  disease 
that  is  characterized  by  a  wide  range  of  symptoms  affecting  many  organ  systems  throughout  the 
body.  Although  larval  zebrafish  have  been  shown  to  often  be  a  good  indicator  of  the 
toxicological  response  in  adult  fish,  we  aimed  to  further  explore  the  effects  PB  and  stress  on 
specific  tissues  by  also  utilizing  adult  fish  for  our  experiments.  We  developed  a  study  to  assess 
both  the  immediate  and  long-term  effects  of  PB  exposure,  with  and  without  co-exposure  to 
stressful  conditions.  To  this  end,  adult  male  zebrafish  (AB  strain,  4-12  months  of  age)  were 
randomly  assigned  to  tanks  to  receive  control,  PB,  stress,  or  PB+stress  treatments.  PB-treated 
fish  were  kept  in  a  5  mmol/L  PB  solution  and  fish  were  stressed  by  30  min  exposure  to  alarm 
pheromone  (35).  A  visual  barrier  was  placed  between  tanks  during  stress  treatment  to  preclude 
alterations  in  behavior  based  on  movement  in  neighboring  tanks.  A  60%  change  of  system 
water  or  PB  solution  was  performed  at  the  conclusion  of  the  30-minute  stress  treatment. 
System  water  served  as  the  control  for  both  PB  and  stress  treatments.  The  fish  in  the 
experiments  were  raised  on  our  recirculating  system,  acclimated  to  static  tanks  for  two  days  and 
received  experimental  treatments  for  three  days.  Following  the  acclimation  and  treatment 
period,  fish  from  half  of  the  tanks  (group  A)  were  euthanized  for  sample  collection  and  the 
remaining  fish  (group  B)  were  used  for  motility  assays  to  assess  neurobehavioral  phenotypes. 
(Figure  13). 

Accftnatinn  WMho^l 
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Figure  13.  Experimental  design  for  study  of  effects  of  PB  and/or  stress  in  adult  zebrafish.  Capital  letters 
indicate  fish  used  for  sample  collection  at  0-wk  (A)  or  8-wk  (B)  post-treatment. 

The  sample  collection  included  excising  a  portion  of  skeletal  muscle  and  the  brain  for  the 
isolation  of  enzymatically-active  protein  and  nucleic  acids  from  each  tissue,  as  well  as  whole 
body  samples  for  determination  of  cortisol  levels.  For  the  motility  analysis,  fish  were  moved 
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individually  from  the  group  treatment  tank  into  100  mL  of  system  water  and  allowed  to  acclimate 
for  15  min  on  the  counter,  followed  by  5  min  of  acclimation  in  the  ZebraBox.  After  the 
acclimation  period,  the  fish  was  subjected  to  the  experimental  protocol  of  light:dark  disruption  as 
described  for  larval  zebrafish.  Fish  used  for  motility  assays  were  moved  to  static  tanks 
containing  system  water  for  a  7-day  washout  period  before  returning  to  the  system.  These  fish 
were  assessed  weekly  for  sustained  changes  in  locomotion  or  behavior,  with  sample  collection 
at  8  weeks  after  the  treatment  period  (Figure  13).  For  the  motility  assay,  this  experimental 
series  was  repeated  twice:  experimental  series  1  consisted  of  four  replications  of  the  protocol, 
experimental  series  2  consisted  of  three  replications  of  the  protocol,  but  with  motility  analysis 
only  at  0,  4  and  8  weeks  post-treatment  to  minimize  the  possibility  of  habituation  being  a 
confounding  factor  in  the  analysis.  In  both  series,  animals  were  between  4  and  12  months  of 
age  and  one  replication  in  each  experimental  series  was  conducted  with  female  fish.  Fish  from 
experimental  series  2  were  euthanized  for  sample  collection  at  12  weeks  post-treatment;  these 
samples  have  been  stored  at  -80°C  and  may  be  used  for  validation  or  follow  up  studies. 

Skeletal  muscle  AChE  activity  was  altered 
immediately  after  treatment  by  an  interaction 
between  PB  and  stress  (p  =  0.058),  with 
stress  attenuating  the  AChE  inhibiting  effect 
of  PB  (Figure  14,  control  vs.  PB  pair  wise  t- 
test  p  =  0.03).  There  were  no  changes 
observed  in  muscle  AChE  activity  at  8-wk 
post-treatment  or  in  brain  AChE  activities  at 
either  time  point.  For  experimental  series  1, 
immediately  following  the  cessation  of 
treatment,  fish  treated  with  PB  traveled  a 
shorter  distance  during  all  phases  of  the 
experiment,  however  a  greater  amount  of 
time  was  spent  in  the  center  of  the  well  and 
specifically  during  the  dark  period  there  were  more  entries  to  the  center  and  a  greater  distance 
traveled  in  the  center  as  well  (Figure  15).  Stress  decreased  the  total  distance  traveled  and 
lowered  the  average  velocity  of  movement  overall  and  in  the  center  of  the  well  during  the  dark 
period,  an  effect  that  was  opposite  of  the  observed  in  larval  zebrafish  and  may  reflect 
differences  in  the  methodologies  used  to  induce  the  stressful  condition  or  an  age-related 
difference  in  the  stress  response.  Replication  of  the  experiment  during  experimental  series  two 
was  designed  to  assess  neurobehavioral  characteristics  only  immediately  after  treatment  and  at 
4  and  8  weeks  post-treatment  to  minimize  the  possibility  of  habituation  occurring  with  increasing 
repetitions  of  the  protocol  and  thereby  confounding  the  results.  RNA-seq  libraries  were  prepared 
and  sequenced  for  brain  and  skeletal  muscle  collected  immediately  after  the  exposure  period  (n 
=  4  per  treatment  group)  and  at  8-wk  post-treatment.  The  increase  in  distance  traveled  in 
center  as  percent  total  at  all  phases  was  also  replicated  in  experimental  series  two.  The 
decrease  in  total  distance  and  increase  in  duration  spent  in  the  center  during  baseline  was  also 
replicated.  Stress  also  resulted  in  a  significantly  decreased  center  velocity  and  there  tended  to 
be  decreased  total  distance  traveled  in  dark  phase  in  fish  from  experimental  series  two.  The 
data  from  experimental  series  one  weeks  one  through  seven  was  analyzed  using  a  linear 
regression  model  taking  into  account  the  main  effects  and  interactions  of  PB,  stress,  and  time 
post-treatment.  In  this  analysis,  there  remained  significant  effects  of  PB  treatment,  stress,  and 
several  instances  in  which  a  PB  x  stress  interaction  was  observed  (Table  2).  PB  treatment 
resulted  in  a  significant  increase  in  the  duration  and  distance  traveled  in  center  as  percent  total 
during  the  first  light  period.  Exposure  to  alarm  pheromone  resulted  in  an  increase  in  the 


control  stress  PB  PB-^stress 


Figure  14.  Acetylcholinesterase  activity  was  inhibited  to  within 
therapeutic  levels  in  the  PB-treated  muscle  at  0-d  post-treatment- 
Data  are  means  {n  =  4)  ±  s.e. 
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Figure  15.  Locomotion  ana  behave  are  atered  n  adult  zeeafish  to'awing  72-h  PB 
and/or  aarm  p/veromone  exposure  during  (A)  base  ne  light  phase.  |B)  dark  stimulus, 
and  (C)  post-stimulus  iignt  ona.se. .  Data  shown  are  from  exoenmenta  senes  t ,  cata 
are  means  (n  =  1 6)  t  s.e. 


Table  2.  An  interaction  between  PB  and  stress  modulates  motility  and  behavior  in  adult  zebrafish  in  the  weeks  after  exposure. 


Total  distance  tmm) 

Center  distance 

(rrm) 

Center, 'total  d  stance 

Center  duration  (s) 

Center  average 

vetocjty  irrnVs) 

Center  entnes 

(COl*p1) 

Baseline  light  phase 

control 

stress 

PB 

PB*  stress 

14977 1 1139*** 
13064  1  1206 
14361  ±899 
11783  ± 1104 

3499 ±  268 
3786  1  474 

36831  344 

38721431 

0.27 1 0.03* 

0.31 1 0.05 

0.31  1  0.04 

0.38  i  0.04 

103.61 10.9* 
114.3H3.9 
108.71 12.1 
121.81 13.4 

39.3 1  3.3”* 

33.6 1 3.0 

36.4  1 2.4 
34.912.9 

94.4  i  5.7 

102.4  ±11.8 
110.319.1 
111.5112.5 

Stimulus  dark  phase 

control 

stress 

PB 

P8*stress 

14646  ±621— 
12622  1694 
130031630 
115101611 

1920  1 120*” 
17261  125 
16591 107 
23871  194 

0.14  10.01” 
0.1510.01 

0.14  i  0.01 

0.25 1  0.04 

45.6  1  5.2”* 
47.614.8 
■14.3 1 5.0 

76.7i  11.6 

46.9  1 2.4”* 

39.6  1  2.4 

42.0  1 2.5 

37.7 1 2.2 

62.0  1  3.5— 

58.1  14.3 

55.8 1  3.8 

76.0 1 4.8 

Post-stimulus  light  phase 

control 

stress 

PB 

PB-«re?s 

13756  1  693*** 
126531  1230 
13211 1867 
126321967 

46181419 
4641  1504 

40161348 

vjLtm 

D.3B  t  0.04 

0.43  t  0.05 

0.34  t  0.04 
P.:39l0-Q4 

139.61  12.3 
143.31 12.8 

122.4  111.1 
133,6-1  HQ 

35.7  1  2.4”* 
34.213.0 
34.612.0 

36.0 All 

116.0  ±10.8 
116.7H2.3 
111.018.9 
ll«-Q.*B.O 

*  p  <  0.05  for  man  effect  of  PB 
”  p  <  0.05  tor  train  effect  of  stress 
p  <  0.05  for  PB  x  stress  interaction 
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percentage  of  distance  traveled  in  the  center  during  the  dark  phase.  In  both  experimental 
series,  in  the  weeks  after  treatment,  in  the  dark  phase  there  were  divergent  behavioral 
responses  to  PB  based  on  the  presence  or  absence  of  stressful  exposure.  PB+stress 
increased  entries  to  center  as  well  as  center  distance  and  duration  and  decreased  total  distance 
and  center  average  velocity,  whereas  the  opposite  effect  was  observed  in  animals  only  treated 
with  PB.  We  have  found  in  these  two  experimental  series  that  PB  and  stress  exposures  alter 
motility  and  behavioral  phenotypes  in  adult  zebrafish;  the  specific  effects  of  PB  treatment 
included  decreased  total  movement  distance  and  duration  immediately  following  treatment, 
whereas  PB  treatment  and/or  PB+stress  treatment  produced  a  more  complex  neurobehavioral 
phenotypes  for  weeks  after  treatment  as  evidenced  by  altered  patterns  in  movement  in  the 
center  of  the  field. 

Libraries  for  mRNA-seq  (n  =  4  per  treatment,  male  fish)  were  sequenced  as  50-bp  paired-end 
reads  on  an  lllumina  HiSeq,  producing  at  least  20  million  reads  for  each  sample.  Reads  were 
aligned  to  the  transcriptome  using  Tophat,  raw  counts  were  determined  using  HTSeq-count, 
and  differential  expression 
of  genes  was  assessed 
using  DESeq  with  a  false 
discovery  rate  (fdr)  cutoff 
of  <  0.05.  In  both  tissues 
and  at  both  0-d  and  8-wk 
timepoints,  we  identified 
sets  of  differentially 
expressed  genes  in 
response  to  PB,  stress, 
and  a  PB  x  stress 
interaction  (Table  3). 

A  targeted  inspection  of  the  data  revealed  that  several  genes  related  to  acetylcholine  signaling 
were  dysregulated.  Subunits  of  the  nicotinic  acetylcholine  receptor  were  downregulated  in  both 
the  brain  and  muscle  immediately  after  the  treatment  period.  We  also  observed  downregulation 
of  choline  acetyltransferase  in  the  muscle  and  of  several  synaptic  proteins  in  the  brain  at  this 
time  point.  Many  of  the  same  or  related  genes  were  differentially  expressed  at  8-wk  post¬ 
treatment,  although  there  was  an  interaction  between  PB  and  stress  at  that  time  point  and  the 
pattern  of  expression  among  treatment  groups  was  not  consistent. 

There  was  minimal  overlap  between  the  specific  genes  in  each  differentially  expressed  gene 
sets,  however,  when  the  DAVID  Functional  Analysis  Tools  were  used  to  identify  functional 
enrichment  there  were  several  instances  of  overlap  in  the  functional  enrichment  categories  for 
PB  differentially  expressed  genes  in  both  tissues  at  the  0-wk  time  points  in  broad  functional 
categories  such  as  regulation  of  transcription,  iron  ion  binding,  and  oxidation  reduction  (Table 
4).  Further  analysis  of  the  genes  found  in  these  categories  may  yield  novel  insights  into  the 
early  systemic  pathogenesis  of  Gulf  War  Illness.  Flowever,  there  were  also  instances  of 
functional  overlaps  at  both  0-d  and  8-wk  post-treatment  within  a  particular  tissue.  In  muscle, 
these  terms  were  related  to  cytoskeletal  structure,  specifically  intermediate  filament,  gap 
junction  and  connexins  for  PB-responsive  genes.  In  brain,  the  extracellular  matrix  proteins  were 
overrepresented  at  both  time  points.  In  addition  to  those  genes  perturbed  by  PB  treatment, 


Muscle 

0-wk  post-treatment 

8-wk  post-treatment 

PB 

306 

303 

stress 

80 

67 

PB  x  stress 

39 

879 

Brain 

0-wk  post-treatment 

8-wk  post-treatment 

PB 

427 

131 

stress 

273 

204 

PB  x  stress 

399 

323 

Table  3.  Differentially  expressed  genes  in  the  brain  or  muscle  of  adult 
zebrafish  after  treatment  with  PB  and/or  alarm  pheromone. 
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Table  4.  Overlapping  terms  for  DE  genes  significant  for  an  effect  of  pyridostigmine  bromide  in  adult  zebrafish  tissues.  Data 
are  p-valueswlth  FDR  rate  correction. 


Cateaixv 

Brain 

Muscle 

Term 

0*wK 

B-wk 

0-wk 

0-wk 

GOTERM  BP  FAT 

GO:C055114-oxidalion  reduction 

49.58 

0.11 

GOTERM  MF  FAT 

G0.0005506~iron  ion  binding 

61.99 

0.08 

SP  PIR  KEYWORDS 

neme 

66.75 

1.28 

SP  PIR  KEYWORDS 

iron 

64.78 

0.29 

GOTERM  MF  FAT 

GO.00037C0-transcripLon  factor  activity 

0.01 

55.62 

SP  PIR  KEYWORDS 

developmental  protein 

0.05 

52.13 

SP  PIR  KEYWORDS 

dna-bindma 

<0.01 

50.82 

GOTERM  CC  FAT 

G0.0005856-cvtoskeleton 

16.58 

3.44 

GOTERM  CC  FAT 

G0.0005882~intormediate  filament 

0.02 

53.77 

GOTERM  CC  FAT 

GO  0005921-gap  lunction 

53.18 

0.36 

GOTERM  CC  FAT 

G0.0005922-connexon  oomglox 

47.50 

0.25 

GOTERM  CC  FAT 

G0.0044430-cytoskeleta  part 

5.33 

4.82 

GOTERM  CC  FAT 

GO.C045 1 1 1  -intermediate  filament  cytoskelelon 

0.02 

53.77 

INTERPRO 

IPR0C0500.Connexjns 

54.41 

0.46 

INTERPRO 

IPR0C1664.hte'med  ate  filament  protein 

0.02 

6B.3B 

INTERPRO 

IPR013092:Connexm.  N-terminal 

52.33 

0.41 

INTERPRO 

IPR01 6044. Filament 

0.02 

63.38 

INTERPRO 

IPR017990:Connexm.  conserved  site 

50.20 

0.36 

INTERPRO 

IPR017991:Connexm  region 

54.41 

0.46 

INTERPRO 

IPR019570.Gap  unction  Drotoin.  cvsteine-ncn  domain 

54.41 

0.46 

SMART 

SM30337:CNX 

29.28 

0.49 

SP  PIR  KEYWORDS 

cel!  (unction 

49.38 

2.48 

SP  PIR  KEYWORDS 

aap  luncion 

45.89 

0.39 

SP  PIR  KEYWORDS 

Intermediate  ‘I'ament 

0.01 

61.79 

GOTERM_CC_FAT 

GO.G005576-extracelular  region 

3.92 

39.80 

Table  5,  Overlapping  terms  for  DE  genes  signif.cant  for  PB  x  stress  interacton 

Data  are  c-values  with  FDR  rate  correction. 
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<0.01 
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IPROOOGSeAnutm 
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SP  PR  KEYWORDS 
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39.27 

0.02 

SP  PR  KEYWORDS 
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11  B3 

1.13 

SP  PR  KEYWORDS 

dna-^nina 

0.01 

<0.01 

SP  PTR  KEYWORDS 

Ho-wbcw 

•0  01 

<0  01 

SP  PR  KEYWORDS 
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3  37 

V  J  ■■ 

DNAT«ningros»a-  Ho^eoPax 
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there  was  enrichment  for  genes  related  to  sensory  perception,  hormone  activity,  and  regulation 
of  transcription  in  expression  altered  by  a  PB  x  stress  interaction  (Table  5). 

Consistent  with  reports  of  crosstalk  and  similar  molecular  perturbation  as  a  result  of  PB  and 
stressful  exposures,  the  steroid  biosynthesis  pathway  was  common  to  several  of  our 
differentially-expressed  gene  lists.  In  the  muscle,  steroid  and  lipid  biosynthetic  processes  were 
downregulated  by  PB  treatment  at  0-d  post-treatment.  Steroid  biosynthesis  was  also 
stronglydownregulated  by  the  stress  treatment  in  both  the  brain  and  muscle  at  the  end  of  the 
72-hr  treatment  period.  However,  a  PB  x  stress  interaction  was  observed  for  this  process  in  the 
muscle  at  8-wk  post-treatment,  with  gene  expression  levels  being  significantly  greater  in 
animals  treated  with  alarm  pheromone  only  compared  to  the  other  groups. 

There  was  enrichment  for  a  number  of  biologically  interesting  term  categories  in  the  PB  x  stress 
interaction  gene  sets  from  muscle  at  8-wk  post-treatment,  such  as  voltage-gated  ion  channel 
activity,  neurotransmitter  transport,  clathrin-coated  vesicle,  and  cell  junction  (PB,  stress,  and  PB 
x  stress  enrichment  terms  found  in  only  one  tissue  and  time  point  found  in  Supporting  Data). 
Enrichment  was  also  observed  for  the  cell  junction  in  the  muscle  at  8-wk  post-treatment  for  the 
genes  that  were  downregulated  in  response  to  PB  treatment.  At  8-wk  post-treatment,  we 
identified  several  functional  categories  that  were  enriched  in  differentially  expressed  genes  from 
both  brain  and  muscle  in  response  to  PB  treatment  or  a  PB  x  stress  interaction,  specifically,  the 
INTERPRO  categories  for  claudin  and  low  density  lipoprotein-receptor  class  A  proteins,  as  well 
as  the  tight  junction  GO  term  and  KEGG  pathway  for  cell  adhesion  molecules. 

These  functional  categories  that  have  been  found  to  be  over-represented  in  our  differentially- 
expressed  gene  sets  have  plausible  biological  connections  to  the  neuromuscular  symptoms  of 
Gulf  War  Illness  and  warrant  further  validation.  These  insights  also  offer  promise  for  the 
development  of  therapies  for  the  treatment  of  Gulf  War  Illness. 


Technical  objective  3:  assessment  of  off-target  effects  by  genetic  and  chemical  means 

To  complement  our  studies  of  PB  in  the  presence  or  absence  of  stressful  conditions  in  AB 
zebrafish  larvae  with  a  genetic  approach,  we  have  also  performed  experiments  to  study  the 
response  to  PB  and/or  stress  in  the  achetm205  and  achetf222a  mutant  strains.  This  line  of  inquiry 
may  suggest  whether  the  effects  of  PB  are  intrinsic  to  changes  in  acetylcholinesterase  activity 
levels  or  whether  they  may  be  attributed  primarily  to  off-target  effects.  Pooled  groups  (n=20)  of 
achetm205  or  achetf222a  larvae  from  clutches  produced  by  either  homozyogous  wild-type  or 
heterzyogous  carriers  were  randomized  to  receive  stress  treatments  with  the  experimental 
protocol  as  described  for  AB  larvae.  The  background  strain  for  both  mutants  is  the  AB  strain, 
however  the  achetf222a  genotype  has  a  point  mutation  adjacent  to  the  catalytic  site  of  AChE, 
whereas  the  achetm205  mutation  creates  a  premature  stop  codon,  resulting  in  a  protein  lacking 
both  the  catalytic  site  and  the  C-terminal  domain  (36).  Both  mutations  were  reported  to  abolish 
AChE  enzymatic  activity.  Given  that  both  of  the  ache  mutations  are  lethal  at  5  dpf,  two-thirds  of 
the  larvae  from  a  cross  of  two  heterozygous  adults  would  be  expected  to  be  two-thirds 
heterozygous  carriers  of  the  mutation.  Assuming  a  complete  lack  of  activity  from  the  mutant 
allele,  no  increase  in  the  expression  of  AChE  and  complete  penetrance,  it  would  be  expected 
that  AChE  activity  would  be  decreased  by  -33%  in  samples  from  pooled  heterozygous  cross 
offspring  samples  and  this  is  what  we  have  observed  in  the  achetm205  strain  (Table  6).  AChE 
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activity  was  significantly  decreased  to  25%  in  the  achetf222a  heterozygous  larvae;  there  was  no 
effect  of  stress  observed  in  any  strain. 


Table  6.  PB  exposure  or  heterozygosity  for  ache  mutant  allele  decreased  acetylcholinesterase 
activity  in  larval  zebrafish.  Data  are  means  ±  s.e. 


Strain 

Inhibition  relative  to  control  (mean  ±  s.e.) 

p-value 

AB 

40.1  ±3.4% 

<  0.01 

ache'”205 

34.0  ±  9.9% 

0.038 

ache™23 

25.4  ±  2.9% 

0.022 

To  further  characterize  which  mutant  strain  produces  a  phenotype  most  similar  to  PB-treated 
AB  larvae,  we  also  conducted  the  motility/neurobehavioral  assay  (as  in  Figure  7)  and 
genotyped  each  larvae  to  identify  heterozygous  carriers.  When  compared  to  the  strongest 
phenotype  observed  in  PB-treated  AB  larvae,  i.e.  decreased  distance  and  duration  in  the  center 
of  the  well  during  the  post-stimulus  period,  we  found  that  there  was  a  significant  decrease  in  the 
distance  traveled  and  a  trend  for  decreased  duration  (p  =  0.084)  in  achetm205'/+  larvae  compared 
to  wild  type  larvae,  whereas  no  differences  were  observed  in  achetf222a/+  larvae  (Figure  16). 
This  is  consistent  with  the  motility  data  from  the  initial  identification  and  characterization  of  the 
abnormal  movement  in  the  ache  mutants  that  found  most  pronounced  phenotype  in  achetm205 
homozygotes  (37).  This  data  suggests  that  the  achetm20S  mutation  produces  phenotypes  that 
are  most  similar  to  the  phenotypes  resulting  from  PB  treatment.  By  comparing  the  chemical 
treatment  to  the  genetic  model  at  the  molecular  level,  we  could  determine  whether  the  effects  of 
PB  are  through  its  primary  mechanism  of  acetylcholinesterase  inhibition  or  whether  off-target 
effects  that  could  contribute  to  the  development  of  Gulf  War  Illness. 


strain  AB 


strain.  aoet'nZSS 

yrslM  -  UWi 


strain  a;neP222a 


Figure  1 6.  AB  larvae  treated  with  PB  arc  heterozygous  ache  --  "  arvae  spent  traveled  less  distance  and/or  spent  less  time  in  the  center  of  the 
well  during  the  post-stimulus  period  in  the  motiity  analysis,  whereas  no  difference  m  behavorwere  noted  ecbe’7*'*.  Data  are  means  {n  =  14-36) 

±  s.e. 


Given  the  possibility  that  off-target  effects  of  PB  may  be  related  to  the  symptoms  of  Gulf  War 
Illness,  it  may  be  desirable  to  consider  other  AChE  inhibitors  (AChEls)  as  potential  alternatives 
for  pretreatment  during  the  threat  of  nerve  agent  exposure.  With  this  in  mind,  we  chose  to  study 
the  response  to  galantamine  (GAL,  brand  name  Razadyne),  rivastigmine  (RIV,  brand  name 
Exelon),  and  donepezil  (brand  name  Aricept).  These  chemicals  are  all  water-soluble  AChEls 
that  are  currently  approved  by  the  FDA  for  the  treatment  of  Alzheimer’s  disease  and  have  been 
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studied  for  their  use  against  nerve  agents.  Larvae  were  treated  with  0.1  mM 
acetylcholinesterase  inhibitor  solution  or 
water  for  72  hrs.  Exposure  to  the 
donepezil  solution  resulted  in  100% 
mortality  of  larvae  within  24-h  after 
administration,  the  experiment  was 
continued  with  the  remaining  treatment 
groups.  At  the  end  of  the  treatment  period, 
protein  and  RNA  samples  were  prepared, 

AChE  activity  was  quantified,  and  mRNA- 
seq  libraries  were  prepared  and 
sequencing  data  analyzed  as  described  for 
the  larval  studies  under  objective  one.  The 
mRNA-seq  libraries  (n  =  3  per  treatment) 
were  sequenced  as  single-end  36-bp  on  an 
lllumina  GAIIx,  producing  -9-30  million 
reads  alignable  to  Zv9  exons  per  sample. 

The  treatment  of  larval  zebrafish  with 
AChEls  was  successful  in  inhibiting  AChE  activity  (Figure  17).  Acetylcholinesterase  was 
inhibited  to  approximately  the  range  of  clinical  relevance  for  GAL  and  PB  treatments. 
Rivastigmine  treatment  produced  60%  inhibition  of  AChE  activity,  which  was  significantly 
different  from  the  control.  Similar  to  the  pattern  for  AChE  activities,  the  greatest  number  of 
differentially  expressed  transcripts  relative  to  control  was  found  in  response  to  RIV  treatment, 
with  both  GAL  and  PB  also  inducing  alterations  in  gene  expression,  but  to  a  lesser  extent 
(Figure  18).  Genes  shared  by  GAL  and  RIV  were  strongly  enriched  for  heat  shock  proteins, 
glutathione  S-transferase,  short-chain  dehydrogenase/reductase,  carbonic  anhydrase, 
scavenger  receptor  activity,  and  extracellular  matrix.  Genes  differentially  expressed  by  PB  and 
GAL  were  enriched  several  terms  we  have  reported  for  the  PB  response  in  our  other 
experiments,  such  as  neurological  system  processes,  synapse,  nucleotide  binding,  and  ion 
transmembrane  transport.  Genes  unique  to  the  PB  response  in  this  data  set  were  enriched  for 
terms  associated  with  sensory  perception,  specifically  retinal  proteins.  These  results  are 

suggestive  that  there  are  few  off-target 
effects  of  PB.  However,  there  were  a 
large  number  of  genes  unique  to  GAL  and 
RIV  exposure  with  significant 
enrichments.  GAL  genes  were 
overrepresented  for  structural  molecules 
of  the  ribosome  and  mitochondria, 
oxidative  phosphorylation  and  cellular 
cation  homeostasis,  whereas  RIV  genes 
were  enriched  for  beta/gamma  crystallins, 
peptidases,  intermediate  filaments,  EF 
hand,  lipocalin,  metabolism  of  xenobiotics 
by  cytochrome  P450,  proteinase 
inhibitiors,  carboxylesterases,  and 
calreticulin/calnexin.  Interestingly,  the 
enrichments  for  the  RIV  unique  DE  gene 
list  share  some  commonalities  with  those 
found  in  the  early  dose  response  study  for 


GAL  PB 


RIV 


Figure  18.  Treatment  of  larval  zebraf  sh  with  acetylcholinesterase  inhibitors  results  in 
differential  gene  express  or.  Numbers  indicate  oount  of  differentially  expressed 
transcripts  (adjusted  p  <  0.05). 
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Figure  17.  Acetylcholinesterase  activity  in  larval  zebrafish  was  inhibited 
by  72-h  treatment  with  acetylcholinesterase  inhibitors  galantamine  (GAL), 
pyridostigmine  bromide  (PB).  and  rivastigmine  (RIV).  Data  are  means  * 
s.e. 
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the  largest  dose  of  PB,  which  produced  levels  of  AChE  inhibition  similar  to  that  produced  by  RIV 
treatment.  In  sum,  the  effects  of  PB  on  gene  expression  appear  to  be  a  secondary 
consequence  of  AChE  inhibition  as  evidenced  by  similar  perturbation  of  gene  expression  in 
response  to  other  AChEls. 


Technical  objective  4:  exposure  offish  to  chemical  with  screening  for  resistance 

In  preparation  for  a  screening  of  the  offspring  of  mutagenized  fish  for  increased  resistance  to 
the  effects  of  PB,  we  exposed  larvae  to  suprapharmacological  PB  doses  (0.01-0.5  mol/L)  to 
determine  the  dosage  at  which  PB  exposure  is  lethal  and  characterized  the  resulting  movement 
phenotypes  at  sublethal  doses.  AB  larvae  were  placed  individually  in  the  wells  of  a  96-well 
plate  and  treated  with  200  pL  of  PB  solution  (0.01,  0.025,  0.05,  0.1,  0.25  or  0.5  mol/L)  with 
randomization  by  column  (n=32  larvae  per  treatment).  After  24  hours  of  treatment,  one  hour  of 
tracking  under  normal  light  conditions  was  used  to  determine  movement  patterns  and  was  used 
in  conjunction  with  visual  assessment  to  determine  the  number  of  larvae  surviving.  While  66% 
of  larvae  survived  at  0.05  mol  PB/L,  only  1  fish  survived  the  0.1  mol  PB/L  dose  and  none  of  the 
fish  exposed  to  either  0.25  or  0.5  mol/L  were  alive  after  24  hours  of  treatment.  As  the  dose  of 
PB  increased,  the  duration  and  distance  of  movement  of  surviving  fish  decreased.  By  tracking 
movement  as  well  as  measuring  the  survival,  it  would  be  possible  to  develop  a  more  sensitive 
screen  in  which  attenuation  of  the  movement  phenotype  could  also  be  utilized  as  an  indicator  of 
resistance  to  PB. 


KEY  RESEARCH  ACCOMPLISHMENTS 

Over  the  duration  of  this  award  we: 

•  Built  an  aquaculture  facility  with  active  breeding  colonies  of  5  strains  of  zebrafish. 

•  Completed  dose-response  and  time-course  studies  of  PB  exposure  in  larval  zebrafish. 

•  Studied  the  effect  of  natural  variation  on  the  response  to  PB  in  larval  zebrafish  by 
comparing  phenotypes  and  dose-responsive  gene  expression  patterns  in  three  strains. 

•  Determined  the  effect  of  stress  on  motility  and  behavioral  phenotypes,  as  well  as  gene 
expression  changes  in  PB-treated  larval  and  adult  zebrafish. 

•  Utilized  ache  mutants  as  a  genetic  model  and  additional  acetylcholinesterase  inhibitors 
as  a  chemical  model  for  studying  the  potential  off-target  effects  of  PB. 

•  Identified  whole-body  and  tissue-specific  perturbations  of  gene  expression  in  response 
to  PB  treatment  using  larval  and  adult  zebrafish  respectively. 


REPORTABLE  OUTCOMES 

A  manuscript  is  in  preparation,  with  submission  expected  by  the  end  of  2013. 
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CONCLUSION 

In  summary,  we  have  utilized  motility/behavioral  phenotyping  and  transcriptional  profiling  in 
larval  and  adult  zebrafish  to  increase  understanding  of  the  effects  of  pyridostigmine  bromide 
and  factors  of  relevance  to  the  combat  situation  in  the  1991  Gulf  War,  such  as  stress  and 
genetic  variation.  Our  data  indicate  that  PB  exposure  alters  behavior  in  larval  and  adult 
zebrafish  as  indicated  by  altered  place  preference,  although  in  larval  zebrafish  the  response 
tended  to  be  anxiogenic,  whereas  an  anxiolytic  effect  was  observed  in  adults.  Our  analysis  of 
the  molecular  response  to  PB,  as  measured  using  mRNA-seq,  repeatedly  identified  genes 
related  to  the  synapse,  calcium  binding  and  signaling,  cytoskeletal  parts,  and  cell  junctions  as 
downregulated.  Immediately  after  treatment  with  the  drug,  we  also  consistently  observed 
perturbations  in  genes  affecting  oxidation-reduction  reactions,  regulation  of  transcription,  lipid 
metabolic  processes,  and  iron  ion  homeostasis.  Similar  categories  of  genes  were  differentially- 
expressed  in  larval  zebrafish  treated  with  the  acetylcholinesterase  inhibitors  galantamine  and 
rivastigmine,  suggesting  that  many  of  these  alterations  of  gene  expression  occur  specifically 
due  to  the  downstream  effects  of  lowered  acetylcholinesterase  activity  and  that  PB  does  not 
have  broad  off-target  effects.  Identification  of  these  categories  of  differentially-expressed  genes 
may  offer  potential  targets  for  therapies  for  the  treatment  Gulf  War  Illness  or  aid  in  the 
prevention  of  the  development  of  similar  pathologies  in  future  conflicts. 
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Functional  annotation  of  DE  genes  for  PB 
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Functional  annotation  of  DE  genes  for  PB 
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Term 
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Functional  annotation  of  DE  genes  for  stress 
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Category 

Term 
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8-wk 

COG  ONTOLOGY 

Secondary  metabolites  biosynthesis,  transport,  and  catabolism 
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G0:0006355~regulation  of  transcription,  DNA-dependent 
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G0:0045449~regulation  of  transcription 
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GO:0003677~DNA  binding 
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one  donor,  and  incorporation  of  one  atom  of  oxygen 

2.58 
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Functional  annotation  of  DE  genes  for  stress 

Category  Term 
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IPR015650:Fleavy  chain  of  Myosin 

<0.01 

SMART 

SM00015:IQ 

0.02 

SMART 

SM00242:MYSc 

0.04 

SP  PIR  KEYWORDS 

actin-binding 

0.04 
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Functional  annotation  of  DE  genes  for  stress 

Brain 

Muscle 

Category 

Term 

0-wk  8-wk 

0-wk 

8-wk 

SP  PIR  KEYWORDS 

coiled  coil 

0.42 

SP  PIR  KEYWORDS 

motor  protein 

0.04 

SP  PIR  KEYWORDS 

muscle  protein 

0.01 

SP  PIR  KEYWORDS 

myosin 

0.09 

SP  PIR  KEYWORDS 

thick  filament 

<0.01 

SMART 

SM00338:BRLZ 

3.30 

Functional  annotation  of  DE  genes  for  PB  x  stress  interaction 

Brain 

Muscle 

Category 

Term 

0-wk  8-wk 

0-wk 

8-wk 

GOTERM  BP  FAT 

G0:0008283~cell  proliferation 

3.09 

INTERPRO 

IPR001699:Transcription  factor,  T-box 

2.19 

INTERPRO 

IPR018186:Transcription  factor,  T-box,  conserved  site 

2.19 

SMART 

SM00425:TBOX 

2.72 

GOTERM  BP  FAT 

G0:0006350~transcription 

3.66 

GOTERM  MF  FAT 

GO:0004703~G-protein  coupled  receptor  kinase  activity 

3.84 

GOTERM  MF  FAT 

G0:0050254~rhodopsin  kinase  activity 

1.58 

INTERPRO 

IPR000239:GPCR  kinase 

3.42 

INTERPRO 

IPR001827:Homeobox  protein,  antennapedia  type,  conserved  site 

<0.01 

PIR  SUPERFAMILY 

PIRSF000584:G  protein-coupled  receptor  kinase 

3.23 

PIR  SUPERFAMILY 

PIRSF002612:homeotic  protein  HoxA5/D4 

<0.01 

PIR  SUPERFAMILY 

PIRSF500608:homeotic  protein  HoxA5 

1.00 

UP  SEQ  FEATURE 

short  seguence  motif:Antp-type  hexapeptide 

<0.01 

INTERPRO 

IPR001 483:Urotensin  II 

3.96 

PIR  SUPERFAMILY 

PIRSF001793:urotensin  II 

3.52 

GOTERM  BP  FAT 

G0:0006694~steroid  biosynthetic  process 

2.52 

GOTERM  BP  FAT 

G0:0006811~ion  transport 

<0.01 

GOTERM  BP  FAT 

G0:0006812~cation  transport 

<0.01 

GOTERM  BP  FAT 

G0:0006813~potassium  ion  transport 

3.51 

GOTERM  BP  FAT 

G0:0006820~anion  transport 

1.51 

GOTERM  BP  FAT 

G0:0006836~neurotransmitter  transport 

3.16 

GOTERM  BP  FAT 

G0:0008202~steroid  metabolic  process 

1.13 

GOTERM  BP  FAT 

G0:000861 0-lipid  biosynthetic  process 

3.33 

GOTERM  BP  FAT 

G0:0015672~monovalent  inorganic  cation  transport 

<0.01 

GOTERM  BP  FAT 

GO:001 5698-inorganic  anion  transport 

3.62 

GOTERM  BP  FAT 

GO:001 6125-sterol  metabolic  process 

0.19 

GOTERM  BP  FAT 

GO:001 6126-sterol  biosynthetic  process 

0.43 

GOTERM  BP  FAT 

G0:0030001~metal  ion  transport 

<0.01 
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Functional  annotation  of  DE  genes  for  PB  x  stress  interaction 

Category  Term 

Brain 

0-wk  8-wk 

Muscle 

0-wk  8-wk 

GOTERM  CC  FAT 

G0:0005887~integral  to  plasma  membrane 

0.27 

GOTERM  CC  FAT 

G0:0008021~synaptic  vesicle 

<0.01 

GOTERM  CC  FAT 

G0:0008076~voltage-gated  potassium  channel  complex 

4.98 

GOTERM  CC  FAT 

G0:0016021~integral  to  membrane 

2.90 

GOTERM  CC  FAT 

G0:0016023~cytoplasmic  membrane-bounded  vesicle 

1.91 

GOTERM  CC  FAT 

G0:0030135~coated  vesicle 

0.38 

GOTERM  CC  FAT 

G0:0030136~clathrin-coated  vesicle 

0.08 

GOTERM  CC  FAT 

G0:0031224~intrinsic  to  membrane 

1.55 

GOTERM  CC  FAT 

G0:0031226~intrinsic  to  plasma  membrane 

0.47 

GOTERM  CC  FAT 

G0:0031410~cytoplasmic  vesicle 

4.68 

GOTERM  CC  FAT 

G0:0031988~membrane-bounded  vesicle 

1.91 

GOTERM  CC  FAT 

G0:0034702~ion  channel  complex 

0.04 

GOTERM  CC  FAT 

G0:0034703~cation  channel  complex 

0.02 

GOTERM  CC  FAT 

G0:0034705~potassium  channel  complex 

4.98 

GOTERM  CC  FAT 

G0:0044456~synapse  part 

<0.01 

GOTERM  CC  FAT 

G0:0045202~synapse 

<0.01 

GOTERM  MF  FAT 

G0:0005216~ion  channel  activity 

<0.01 

GOTERM  MF  FAT 

G0:0005230~extracellular  ligand-gated  ion  channel  activity 

2.60 

GOTERM  MF  FAT 

G0:0005244~voltage-gated  ion  channel  activity 

0.03 

GOTERM  MF  FAT 

G0:0005261~cation  channel  activity 

<0.01 

GOTERM  MF  FAT 

GO:001 5267-channel  activity 

<0.01 

GOTERM  MF  FAT 

G0:0015293~symporter  activity 

0.47 

GOTERM  MF  FAT 

G0:0022803~passive  transmembrane  transporter  activity 

<0.01 

GOTERM_MF_FAT 

G0:0022832~voltage-gated  channel  activity 

0.04 

GOTERM  MF  FAT 

G0:0022836~gated  channel  activity 

<0.01 

GOTERM  MF  FAT 

G0:0022838~substrate  specific  channel  activity 

<0.01 

GOTERM  MF  FAT 

G0:0022843~voltage-gated  cation  channel  activity 

0.06 

GOTERM  MF  FAT 

G0:0030955~potassium  ion  binding 

0.23 

GOTERM  MF  FAT 

G0:0031420~alkali  metal  ion  binding 

0.03 

GOTERM  MF  FAT 

G0:0046873~metal  ion  transmembrane  transporter  activity 

<0.01 

INTERPRO 

IPR000532:Glucagon/GIP/secretinA/IP 

1.00 

INTERPRO 

IPR001664:lntermediate  filament  protein 

<0.01 

INTERPRO 

IPR003091  :Voltage-dependent  potassium  channel 

3.68 

INTERPRO 

IPR004031  :PMP-22/EMP/MP20/Claudin 

1.99 

INTERPRO 

IPR005821  don  transport 

0.39 

INTERPRO 

IPR006821  dntermediate  filament,  DNA-binding  region 

0.09 
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Functional  annotation  of  DE  genes  for  PB  x  stress  interaction 

Category  Term 

Brain 

0-wk  8-wk 

Muscle 

0-wk  8-wk 

INTERPRO 

IPR017995:Homeobox  protein,  antennapedia  type 

<0.01 

KEGG  PATHWAY 

dre00100:Steroid  biosynthesis 

<0.01 

SMART 

SM00070:GLUCA 

0.55 

SP  PIR  KEYWORDS 

cell  junction 

0.06 

SP  PIR  KEYWORDS 

duplication 

4.66 

SP  PIR  KEYWORDS 

ion  transport 

<0.01 

SP  PIR  KEYWORDS 

ionic  channel 

<0.01 

SP  PIR  KEYWORDS 

iron 

3.39 

SP  PIR  KEYWORDS 

lipid  synthesis 

0.09 

SP  PIR  KEYWORDS 

membrane 

0.22 

SP  PIR  KEYWORDS 

potassium 

0.76 

SP  PIR  KEYWORDS 

potassium  transport 

1.12 

SP  PIR  KEYWORDS 

Sodium  transport 

0.04 

SP  PIR  KEYWORDS 

Steroid  biosynthesis 

0.23 

SP  PIR  KEYWORDS 

sterol  biosynthesis 

0.07 

SP  PIR  KEYWORDS 

Symport 

2.81 

SP  PIR  KEYWORDS 

synapse 

0.50 

SP  PIR  KEYWORDS 

Transcription 

0.89 

SP  PIR  KEYWORDS 

transcription  regulation 

1.01 

SP  PIR  KEYWORDS 

transmembrane 

0.01 

SP  PIR  KEYWORDS 

transport 

1.02 

SP  PIR  KEYWORDS 

voltage-gated  channel 

3.66 
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